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Abstract 

The work by Soares et al. [Phys. Rev. B 65, 174506 (2002)] investigates the 
BCS-BE crossover for <i-wave pairing in the 2-dimensional attractive Hub- 
bard model. Contrary to their claims, we found that a non-pairing region 
does not exist in the density vs coupling phase diagram. The gap parameter 
at T = 0, as obtained by solving analytically as well as numerically the BCS 
equations, is in fact finite for any non-zero density and coupling, even in the 
weak-coupling regime. 

PACS numbers: 74.20.-Z, 74.25.-q, 74.25. Jb 

In Ref. IJ, Soares et al. analyzed the BCS-BE crossover at T = for o?-wave pairing in 
the attractive Hubbard model, extending the work by den Hertog by including the next- 
nearest neighbor hopping t' in the tight-binding dispersion. As is evident from the insets of 
Figs. 1 and 3 as well as from the phase diagram of Figs. 2 and 4 of Ref. |IJ, in the weak- 
coupling regime and for any density, they obtained a non-pairing metallic phase characterized 
by an exactly vanishing value of the BCS gap parameter at T = (A(T = 0) = 0). Similar 
results have been previously obtained by den Hertog in Ref. for t' = 0. In a related 
context, a metallic phase at T = in the weak-coupling regime could only result from the 
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Capone et al. calculation in Ref. || (using a DMFT approach to the s-wave attractive 
Hubbard model at T = 0) by suppressing superconducting pairing. On the other hand, 
results for the gap parameter obtained with Hubbard-like lattice models can be connected 
(in the weak-coupling and low-density limits) to the analytic results of the continuous models 
(see below) obtained in Refs. || and 0, where a finite value for the gap is found in the 
weak-coupling regime. 

In this Comment, we show that the (i-wave gap parameter at T = is always finite for any 
non-zero density and even in the weak-coupling regime, by performing both numerical and 
analytic calculations. This result was reported in Ref. where it was explicitly contrasted 
with the results by den Hertog 0; further comments on this point can be found in the 
review article by Loktev et al. 0. In particular, Soares et al. seem not to have been aware 
of the results by Ref. || concerning the weak-coupling regime. The difference in the results 
must be due to inaccuracy in the numerical calculations of Refs. [p]] and [Q] . 

To further support our previous finding about the existence of a (i-wave BCS supercon- 
ducting ground state (A(T = 0) ^ 0) in the weak-coupling regime, we solve the coupled 
self-consistent equations for the gap function and density at T = 0, as given by BCS theory: 



Here, £ k — — 2t(cos k x + cos k y ) — 4i' cos k x cos k y — \i is the tight-binding dispersion with 
nearest- and next-nearest-neighbor hopping, \i the chemical potential, N the total number 
of sites, and Vk,k' — — ^7k7k' (with 7 k = cos k x — cos k y ) the pairing interaction in the d- 
wave channel. (We set the lattice spacing to unity.) The wave- vector dependence of the gap 
parameter resulting from Eq. (1) is purely (i-wave, i.e., A k = A7 k , provided the interaction 
Vk,k' can be written in a separable form. This considerably simplifies the solution of the 
above equations. 
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We first solve the coupled Eqs.(l) and (2) by the standard Newton's method with high 
precision (using a fc-space mesh up to 1024 x 1024 points), with the value t'/t = —0.1 || to 
compare with the inset of Fig. 3 by Soares et al. Our numerical results for the maximum 
value of the gap parameter A(T = 0) are shown in Fig. 1 as a function of the interaction 
strength V/At and for different densities. In Fig. 1, the arrows locate the values of V/4t 
below which the gap would vanish (for given density) according to the numerical results by 
Soares et al. We find instead finite values for the gap even well below the critical couplings 
identified by Soares et al. (and also by den Hertog || for the case t' = 0). In the weak- 
coupling regime, our result is that the gap parameter decreases exponentially as the coupling 
is decreased, in agreement with the standard BCS result. We obtain the same result below 
by analytic calculations. Regarding instead the chemical potential, our results agree with 
those reported in Fig. 3 by Soares et al., as anyway expected by the smallness of the ratio 
A//i in the weak-coupling region where our results for A disagree from those by Soares et 
al. 

In the low-density regime, where the dispersion is parabolic, the c?-wave BCS equations 
can be solved analytically in the weak-coupling limit. The expansion of the dispersion and 
of the (i-wave factors for small wave vectors leads to £k — {t + 2t')(k 2 , + k 2 ) — /i — 4t — At' = 
k 2 /(2m*)—eF and 7k ~ (k 2 . — k' 2 .)/2 = (k 2 /2) cos 20, respectively, where the polar coordinates 
(A;, 4>) have been introduced together with effective mass m* = 1/(2 (t + 2t')) and Fermi level 
ep = fi+4t+4t'. In addition, in the weak-coupling regime it is convenient to limit the integral 
over the energy variable e to a small window | e — tp \< Uq/2 about the Fermi level, where 
u is a cutoff such that uj ^> A. The integral outside this energy window gives sub-leading 
contributions that can only be evaluated numerically since the full wave-vector dependence 
of £k and 7k ought to be retained in this case. Note that, in the weak-coupling regime, the 
Fermi level coincides with the non-interacting value ep = k F /(2m*) = irn/m* and Eqs. (1) 
and (2) can be decoupled. After these manipulations, we are led to the following expression 
for the <i-wave gap equation at T = 0: 
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7T y-wo/2 Jo J^ + A 2 d COS 2 9 



where g& = Vk%m* / '(167r), A d = Ak 2 F /2, 9 = 2(f), and £ = e — 6p. Equation (3) has been 
written in a form which maps exactly into Eq. (A2) of Ref. ||. The analytic solution for 
the gap equation in the weak-coupling limit reported in Ref. clearly shows that the gap is 
always finite, no matter how weak the effective attraction is. Since this point is crucial to 
the present discussion, we provide here some details of the derivation of the analytic result 
given by Eq. (A2) of Ref. |J. The integral over the angular variable in Eq. (3) can be 
performed exactly using known results for the elliptic integrals 0. One obtains: 

I ,„ cos 2 9 1 rf cos 2 ^ 

ad , = , / ad 



£ 2 + A 2 cos 2 # \/e + A 2 Vl-/« 2 sin 2 # 

(4) 



e + A 2 
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where k 2 = A^/(£ 2 + A 2 ,), and k) and i£(^,Ac) are elliptic integrals of the first and 
second kind, respectively. Performing at this point the integral over the variable £, the 
expansion of the elliptic functions J5(§,«) and up to second-order in k produces 

a log contribution in u /A d , leading to the usual BCS exponential dependence of the gap 
parameter on coupling. The prefactor of this exponential is obtained by considering the sum 
of the remaining terms of the above expansion, resulting in finite integrals over the variable 
£ as ojq/ Ad — > oo. These terms can be summed exactly by exploiting the analytic result 
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The final expression for the gap parameter A^ obtained in this way is: 

A d = ^ exp(-l/^) . (6) 

Note that the coupling constant gd = 7rn 2 Vm* /4 entering the exponential in Eq. (6) depends 
on the squared density, leading to a marked exponential suppression of the gap parameter in 
the weak-coupling regime as the low-density limit (n — > 0) is approached. We have verified 
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that our numerical results for the gap parameter reproduce this dependence on coupling 
strength V and density n. We emphasize, however, that only the exponential dependence in 
Eq. (6) can be compared with the numerical results, since the sub-leading contributions to 
the energy integral in Eq. (3) outside the Uo-Wmdcw (which are not included in the present 
calculation) are expected to eventually modify the prefactor of Eq. (6) in such a way to 
eliminate this cutoff completely (the energy scale ujq is absent in our starting model and has 
been introduced only for convenience in the analytic calculation). 

Summarizing, we have found that a non-pairing region is completely absent in the 
phase diagram of the attractive Hubbard model for c?-wave pairing describing the BCS- 
BE crossover. The gap parameter at T = is, in fact, finite for any non-zero density and 
coupling, even in the weak-coupling regime, contrary to the results of Soares et al. |J and 
den Hertog @. 
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FIGURES 
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FIG. 1. Maximum value of the ci-wave gap parameter A at T = as function of the interaction 
strength V for different densities [n = 0.1 (plus), n = 0.2 (cross), n = 0.3 (star)]. Energies are 
normalized with respect to the half-bandwidth At. The arrows locate the values of V/4t below 
which Ref.fl] claims that A = 0. 
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